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ABSTRACT 

We use Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5) u, g, r, i, z photometry to study Milky 
Way halo substructure in the area around the North Galactic Cap. A simple color cut {g — r < 0.4) 
reveals the tidal stream of the Sagittarius dwarf spheroidal, as well as a number of other stellar 
structures in the field. Two branches (A and B) of the Sagittarius stream are clearly visible in an 
RGB-composite image created from 3 magnitude slices, and there is also evidence for a still more 
distant wrap behind the A branch. A comparison of these data with numerical models suggests that 
the shape of the Galactic dark halo is close to spherical. 

Subject headings: galaxies: kinematics and dynamics — galaxies: structure — Local Group — 
Sagittarius dSph - Milky Way: halo 



1. INTRODUCTION 

Stellar streams in the Milky Way halo produced 
by the accretion of smaller galaxies are a stan- 
dard pred iction of hierarchical mer ging cosmogo- 
nies (e.g., iLvnden-Bell fc Lvnden-Belil Il99 A and ref- 
erences therein). The most spectacular example is the 
disrupting Sagittarius dwarf spheroidal (Sgr dSph), orig- 
inally discovered by Ibata, Gilmore & Irwin (1995). It 
has a heliocentric distance of ~ 25 kpc and is centered 
at Galactic coord inates of £ = 5.6° and b — —14.0° 
l)Ibata et al.lll99'^ . It is dominated by an intermedi- 
ate age population (between 6 and 9 Gyrs, Bellazzini 
et al. 2006), but there is evidence for a much older pop- 
ulation (> 10 Gyrs) as well ( Monaco ct al. 2003). The 
metallicity [Fe/H] ranges from very metal-poor (as low as 
-2 based on the globular clusters) up to approximately 

solar , with probably a mean of 0.5 (|M onaco et aLl 

120051 and references therein). It was realised early 
on that there was some tidal debris in the neighbour- 
hood of the Sgr dSph ( Ibata et al . 1997; Maiewski et al.l 
119991) and that the distribution of this material trac ed the 
Sgr dSph orbit. Subsequently, lYannv et all 1)200(1) used 
Sloan Digital Sky Survey (SDSS) first-year commission- 
ing data to identify an over-density of blue A-type stars 
in two strips located at (^, 6, i?) = (350°, 50°, 46 kpc) 
and (157°, —58°, 33 kpc), which were then matched with 
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the Sgr stream ijlbata et al.ll200lD . Likewise. Ilvezic et alJ 
( 2000) noticed that clumps of RR Lyrae stars in SDSS 
commissioning data lay along the Sgr stream's orbit. 

The be st panorama of the Sgr stream to date was ob- 
tained bv lMaiewski et ^ (2003) using M giants selected 
from the Two Micron All-Sky Survey (2MASS). They 
saw the trailing tidal tail very clearly in the Southern 
Galactic hemisphere, as well as part of the leading arm 
reaching towards the North Galactic Cap. Here, we use 
SDSS Data Release 5 (DR5) to provide a picture of the 
leading arm of the Sgr stream in the vicinity of the North 
Galactic Cap with remarkable clarity, together with a 
number of other notable stellar structures in the field. 

2. THE DATA AND A SIMPLE COLOR CUT 

SDSS llYork et alJ I2OO0I) is an imaging and spec- 
troscopic survey that has mapped ^ 1/4 of the 
sky. Imaging data are produced simultaneously 
in five photometric bands, namely m, q , r, i, and 
z (iFukugita et alJ IT991 iGunn et all 119981: iHogg et al 
2001 : lAdelman-McCarthv et all l2006t iGunn et al 
2000). The data are processed through pipelines 
to measure photometric and ast r ometric proper 



ties (iLupton. Gunn. fc Szalavl Il999t iStoughton et al 



2(tn i iSmith et all l2002t iPier et aljl2003l: Ilvezic et al 
20041) and to select targets for spectroscopic 
follow-up. For de-reddenin g, we use the maps of 
iSchlegel. Finkbeiner. & DavisI l|1998D . DR5 covers 
~ 8000 square degrees around the Galactic North 
Pole, together with 3 strips in the Galactic southern 
hemisphere. We use the catalogue of objects classified as 
stars with artifacts removed together with magnitude 
limits r < 22 and g < 23. At low right ascension and 
declination, we are limited by the boundary of DR5. We 
choose to cut our sample at a = 230° and 6 — 60°. This 
gives us a total of « 2 x 10^ stars. 

Figure n shows the density of stars satisfying the color 
cut g - r < 0.4 in SDSS DR5. As we show shortly, 
these are upper main sequence and turn-off stars belong- 
ing to the Sgr stream. This is a RGB-composite im- 

See |http://cas.sdss.org/astro/en/help/docs/realquery.asp#flags| 
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Fig. 1. — The spatial density of SDSS stars with g — r < 0.4 around the North Galactic Cap in equatorial coordinates, binned 
0.5 X 0.5 arcdegrees. The color plot is an RGB composite with blue for the most nearby stars with 20.0 < r < 20.66, green for stars with 
20.66 < r < 21.33 and red for the most distant stars with 21.33 < r < 22.0. Note the bifurcation in the stream starting at a fx 180°. 
Further structure that is visible includes the Monoceros Ring a,t a ^ 120° , and a new thin stream at 150° J; « ^ 160° and 0° < <5 < 30° . 
The color bar shows a palette of 50 representative colors labeled according to the stellar density (in units of 100 stars per square degree) 
in each of red, green and blue components. The displayed density ranges are 102 to 330 (red), 107 to 304 (green) and 98 to 267 (blue). 




age assembled from three grayscale images colored red, 
green and blue corresponding to the density of the se- 
lected stars in the following three magnitude bins: red 
is 21.33 < r < 22.0, green is 20.66 < r < 21.33, and 
blue is 20.0 < r < 20.66. For a given stellar popula- 
tion, these magnitude bins are distance bins, with red 
being the most distant and blue the nearest. There is a 
clear distance gradient along the stream, from the nearer 
parts at right ascensions a « 120° to the more distant 
parts at a « 210°. There is also a bifurcation visible in 
the stream starting at about the North Galactic Pole; in 
what follows, we will refer to the lower declination branch 
as branch A, and the higher declination as branch B. 

Majewski et al. (2003) traced the Northern stream of 
the Sgr for right ascensions a between 270° and 190°. 
For a < 190°, Majewski et al. (2003) did not see a clear 
continuation of the stream. The combination of Majew- 



ski et al.'s (2003) M giants together with the SDSS stars 
in Figure shows for the first time the entirety of the 
stream, including its continuation through the Galactic 
Cap and into the Galactic Plane. Figure El also shows 
the locations of a number of globular clusters, some of 
which are known to be associated with the Sgr stream. 
For example, Bellazzini et al. (2003) used 2MASS data 
to conclude that NGC 4147 was physically immersed in 
the stream. 

Figure n displays such a remarkable wealth of Galactic 
substructure that it might appropriately be called the 
"Field of Streams". Among the most visible of these 
is the whitish-blue colored, and hence relatively nearby, 
stellar overdensity centered at (a « 185°, (5 « 0°), found 
by Juric et al. (2005) and named the Virgo Overden- 
sity; this is perhaps the same structure as the nearby 
2MASS Northern Fluff (Majewski et al. 2003). Parts 
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TABLE 1 

Locations of the on-stream fields along branch A (a, 5 a) and branch B (a, (5s). The companion off-stream fields for 

BRANCH A HAVE THE SAME RIGHT ASCENSION BUT ARE OFFSET IN DECLINATION BY +20° FOR a > 150° AND BY —12° OTHERWISE. FOR 

BRANCH B, THE OFF-STREAM FIELDS ARE OFFSET BY +13° IN DECLINATION. 
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Fig. 3. — An RGB composite of the SDSS stars as in Figure 1, but now in Galactic coordinates (£, b). Note in addition to the branches 
of the Sgr stream, a second orphan stream is clearly visible at latitudes of b « 50° and with longitudes satisfying 180° ^ ^ ^ 230° . The 
Monoceros Ring is also discernible at low latitudes. 



of the Monoceros Ring (Newberg at al. 2002) are visi- 
ble as the blue-colored structure at a « 120°. Figure 13 
an RGB-composite image of the SDSS stars in Galac- 
tic coordinates 6), also shows the arc-like structures 
of the Monoceros Ring, as predicted by the simulations 
of Penarrubia et al. (2005). Two of the globular clus- 
ters with tidal tails previously identified in SDSS data 
- namely Pal 5 (Odenkirchen et al. 2001) and NGC 
5466 (Belokurov et al. 2006) - can be discerned in 
the figures, together with their streams. Finally, a new 
stream is shown clearly, running from a,5 k, 160°, 0° to 
a,(5 « 140°, 50° ib « 50° and 180° < I < 230° in Fig- 
ure|3Jl- It is distinct from the Sgr stream, which it crosses; 
we discuss its progenitor in a future contribution. 

3. TOMOGRAPHY OF THE SGR STREAM 

To analyze the three-dimensional structure of the 
stream, we set up a series of 6° x 6° fields along branches 
A and B, shown as red (for A) and black (for B) squares 
in Figure 13 The first three A fields actually probe both 
the A and B branches, which are merging at these loca- 
tions. The coordinates of the field centers are listed in 
Table 1. For each on-stream field, there is a companion 
off-stream field of size 15° x 15°, which has the same right 
ascension but is offset in declination as noted in Table 1. 
The off-stream fields are larger, so that the background 
is as smooth as possible. 

Color-magnitude diagrams (CMDs, g — i versus i) are 
constructed for each of the fields and then normalised 
by the number of stars. The difference between each on- 
stream field and its companion off-stream field reveals 
the population of stars belonging to the Sgr stream. An 
example for the field A7 is shown in Figure^fa), together 
with a one-dimensional slice at g—i = 0.55 in Figure^fb). 
The subgiant branch is clearly visible and its location can 
be found to good accuracy. In fact, we fit a Gaussian to 
the one-dimensional slice to obtain its i band magnitude 



and uncertainty. For the example in Figure 0fa), there 
are two subgiant branches visible in the CMD, and two 
distinct peaks in the luminosity functions in Figure 2fb). 
These correspond to two distinct structures at different 
distances, possibly different wraps of the Sgr stream. Al- 
though they could correspond to different populations at 
the same distance, this seems unlikely as the magnitude 
difference between the two subgiant branches changes on 
moving along the stream. 

The Sgr dSph is known to contain a variety of stellar 
populations with different ages and mctallicities. The 
most recent study of this is bv iBelokurov et al.l (|2006), 
who presented a comprehensive CMD of the body of 
Sgr. Rather than attempting to fit multiple theoretical 
isochroncs to our CMDs, we show a direct comparison 
between Belokurov et al. ( 2006) 's CMD and a composite 
CMD for the whole of the stream. The composite CMD is 
created by using the subgiant branch location to measure 
the magnitude offset. We use fields A4 to A16, with A16 
as the reference. The result, transformed to B — V versus 
V (Smith et al. 2002), is shown as the black contours in 
Figure^l^c). This is overlaid upon the colored contours of 
Bellazzini et al.'s (2006) CMD, corrected for extinction 
in i? — V" and V by -0.12 and -0.4 magnitudes, respec- 
tively. The upper main sequence, turn-off, and subgiant 
branches are all well-matched. The only significant dif- 
ference occurs in the sparsely populated upper red giant 
branch (the dark blue contours). Hence, this stream is 
entirely consistent with being composed of the same mix 
of populations as the Sgr dSph. The 0.6 magnitude offset 
between our composite CMD and Bellazzini et al.'s is a 
direct measure of distance. 

Figure El^d) shows the i band magnitude of the sub- 
giant branch versus right ascension of the fields. The 
triangles and stars correspond to the two structures in 
the A branch, the diamonds to the B branch. Only one 
sequence is detected in the B fields. Although its dis- 
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Fig. 4. — (a): Color-magnitude histogram (Hess diagram) of the on-stream field A7 (that is, the 7th field on the A branch) minus the 
number-scaled corresponding off-stream field. The color bin used to construct the one-dimensional slice is marked, (b): The luminosity 
functions of the one-dimensional slices for the on-stream field (unbroken), the off-stream field (broken) and the difference (dotted). Notice 
the two subgiant peaks corresponding to two structures at different distances or two distinct stellar populations, (c); A composite CMD 
for branch A referenced to field A16 (black contours) overplotted on the CMD from Belokurov et al. (2006) (colored contours) of the body 
of the Sgr dSph. (d): Magnitude versus right ascension of the subgiant branch for the A fields (stars and triangles) and for the B fields 
(squares). The distance in kpc is given on the right-hand side axis. If a datapoint is missing, a detection was not possible in that field. 



tance offset is significant only at about the 1 a level, the 
B branch is systematically brighter and hence probably 
slightly closer. Assuming the distance to the Sgr dSph is 
~ 25 kpc, then Figure^d) can be converted directly into 
heliocentric distance versus right ascension, as shown on 
the right-hand axis. The two distinct structures in the 
A fields are separated by distances up to « 15 kpc. Note 
that there is no evidence for any part of the Sgr stream 
passing close to the Solar neighbourhood, as has some- 
times been conjectured (Frecsc et al. 2004). 

This structure of three branches - two close together 
and one more distant - is understandable on comparison 
with the numerical simulations. For example, the effect 
is clearly visible in the lower panel of Figure 3 of Helmi 
(2004), where there are two close branches representing 
material stripped between 3 Gyr and 6.5 Gyr ago, and 
less than 3 Gyr ago. Branches A and B are therefore 
probably tidal debris, torn off at different times. The 
older material, stripped off more than 6.5 Gyr ago, lies 
well behind the two close branches in Helmi's simula- 
tions. This earlier wrap of the orbit may correspond to 
the distant structure seen behind the A branch. 

4. CONCLUSIONS 

We used a simple color cut g — r < 0.4 to map out 
the distribution of stars in SDSS DR5. The "Field of 
Streams" is an RGB-composite image composed of mag- 
nitude slices of the stellar density of these stars. It reveals 
a super-abundance of Galactic substructure, including 
the leading arm of the Sgr stream, as well as a number 
of sibling streams (some hitherto unknown). 

Part of the Sgr stream has previously been seen in the 
Northern hemisphere by Majewski et al. (2003). Here, 
we have mapped out the continuation of the stream, as 
it passes by the Galactic Cap and returns to the Galactic 
Plane. At least two branches of the stream - labelled A 
and B - have been identified, corresponding to material 
torn off at different epochs. There is also evidence for a 
still more distant structure behind the A branch. 

The Sgr stream provides a probe of the shape of the 
Galactic halo. Previous work has been hampered by the 
absence of data i n the most im portant regions of the 
sky. For example, iHelmi I l)2004(l provides simulations of 



the Sgr stream for a range of halo shapes from extreme 
oblate to prolate, all of which broadly agree with the data 
available at that time. Our data covers the critical re- 
gion where the models differ substantially. Examination 
of the simulations displayed in Figure 2 of Helmi ( 2004|) 
enables some preliminary conclusions to be drawn. In 
oblate haloes, the Sgr stream is much fatter and the stars 
much more scattered than shown in our Figure 1. Only 
the spherical and very mildly prolate haloes (1 < q < 1.1, 
where q is the axis ratio in the potential) seem reason- 
able matches to the data. In both these simulations, the 
Sgr stream shows the same bifurcation and overall mor- 
phology as in the SDSS data. We suggest that this may 
be a powerful discriminant of halo shape (Fellhauer et al. 
2006). 
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